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Abstract. We report the discovery of a substellar companion to the intermediate-mass star HD 11977 (G5III). 
Radial velocities of this star have been monitored for five years with FEROS at the 1.52-m ESO and later at the 
2.2-m MPG/ESO telescope in La Silla, Chile. Based on the collected data we calculated an orbital solution with 
, a period of P = 711 days, a semi-amplitude of K\ = 105ms _1 , and an eccentricity of e = 0.4. The period of the 

■ radial-velocity variation is longer than that of the estimated stellar rotation, rendering it unlikely that rotational 

modulation is the source of the variation in the radial velocity. This hypothesis is supported by the absence of 
, a correlation between stellar activity indicators and radial-velocity variation. By determining a primary stellar 

mass of Mi, = 1.91 Mq, the best-fit minimum mass of the companion and semi-major axis of the orbit are 
m.2 sin i = 6.54 Mj up and a2 = 1.93AU, respectively. An upper limit for the mass of the companion of 
ni2 Ss 65.5 Mjup has been calculated from HlPPARCOS astrometric measurements. Although the possibility of a 
brown-dwarf companion cannot be excluded, HD 11977 B is one of the few planet candidates detected around 
an intermediate-mass star. The progenitor main-sequence star of HD 11977 is probably an A-type star. This 
discovery gives an indirect evidence for planetary companions around early type main-sequence stars. 
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' 1- Introduction Intermediate-mass main-sequence stars have spectral 

classes from A to early F and higher effective tcm- 

Among the substellar companions discovered around perature s than solar-type stars. Therefore, their spectra 

other stars, there are only a few planetary and show fewer absorption lines than later, solar-like spectral 

brown-dwarf companions to intermediate-mass stars typeg (late Fj q and K) Due to faster stellar rota tion 

(M* « 1.5-7 M ). In fact, the presence of planetary or (wsin ■ > aokms' 1 ) their line profiles are also broader, 

brown-dwarf companions to stars in this mass regime These facts make it difficult to determine accurately the 

is not well explored in current radial-velocity planet stellar radial velocity (RV) 
search programmes. Ongoing radial-velocity surveys, e.g. 

CORALIE & HARPS survey and California & Carnegie Cool evolved stars, on the other hand, are suitable 

Planet Search, concentrate on solar-like stars. Only few targets for precise RV measurements. As a star moves 

planetary companions have been detected around evolved towards the red giant branch (RGB), its rotational ve- 

stars, like G and K giants (Hatzes et al. I2UU51 and refer- lo chy decreases and the stellar atmosphere cools signifi- 

ences herein). cantly. Red giants, however, are chromospherically active 

(Pasquini et al. I2UUU|I . The variation of the chromosphcric 

Send offprint requests to: J. Setiawan, setiawan@mpia.de activity of G-K giants leaves an imprint in the spectral 

* Based on observations collected at the 1.52-m ESO tele- line profiles. Rotational modulation due to starspots can 

scope from October 1999 to February 2002 and at the 2.2-m also produce RV variations, thereby mimicking the grav- 

MPG/ESO telescope from December 2003 to November 2004 itational influence of low-mass companions. This effect 

at the La Silla Observatory, Chile. can be investigated through the variations in the spec- 
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tral line profile asymmetry and by monitoring the chro- 
mospheric activity indicators, like Call K and H emission 
lines, see e.g., Choi et al (1995J. Non-radial oscillations 
may also induce RV variations as shown by the astero- 
seismology of several giant stars (Hatzes & Cochran ll994l 
Buzasi et al. lSUUUl Frandsen et al. l2UU2|l . 

In this letter we report the discovery of a substellar 
companion around HD 11977. We have excluded rota- 
tional modulation as the source of the RV variations, as 
described in the following sections. 

2. About the star HD 11977 

HD 11977 is a cool evolved star of spectral type G5III 
(HiPPARCOS, ESA 1997). In the Hertzsprung-Russell di- 
agram (My vs. T c ff) it is located in the "clump region" 
of the RGB. After ascending the RGB, the star moved 
down to the clump region where it is undergoing helium- 
core and hydrogen-shell burning. Clump giants typically 
spend 10 s years in this evolutionary stage, that, for initial 
masses higher than 1.7 Mq, represents more than 5% of 
their main-sequence lifetime (Girardi et al.[2000). 



Table 1. Stellar parameters of HD 11977 



Spectral type 




G5III 




mv 




4.68 


mag 


M v 




0.57 ±0.07 


mag 


B-V 




0.931 ± 0.036 


mag 


Parallax 




15.04 ± 0.47 


mas 


Distance 




66.5 ± 2.1 


pc 


Effective temperature T e ff 


4970 ± 70 


K 


Stellar radius R* 




10.09 ±0.32 


R 


Angular diameter 




1.49 ±0.02 


mas 


Metallicity [Fe/H] 




-0.21 ±0.1 


dex 


Surface gravity log ( 


ff/cms 2 ) 


2.90 ±0.2 




Stellar mass M* 




1.91 ±0.21 


M 


log (Age/year) 




8.7-9.3 


t>rot sin i 




2.4 ± 1.0 


kms 


Micro turbulence 




1.40 - 1.60 


kms 


Prot / sin i 




230 - 270 


days 



The basic stellar parameters of HD 11977 are listed 
in Table 1. The spectral type, photometric data and 
parallax were taken from the HiPPARCOS catalogues. 
Adelman 1|2UU1[1 reported a photometric variability of 
0.5 millimag. Thus, HD 11977 is one of the least variable 
stars known from the HiPPARCOS photometric survey. 

The solar-metallicity evolutionary tracks by 
Girardi et al. l|27)07)| show that the mass of HD 11977 
is, indicatively, between 2.5-3.0 Mq. In a forthcoming 
paper (da Silva et al., in preparation) the T e fj, i?*, [Fe/H] 
and \ogg are measured from the spectra of HD 11977. 
Their values are given in Table 1. Then, using as a 
constraint the observed My, B — V, [Fe/H] values and 
their uncertainties, and the Girardi et al. theoretical 
isochrones, we derived the mass probability distribution 
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Fig. 1. Radial velocity measurements of HD 11977, taken 
with FEROS from October 1999 to November 2004. The 
measurements can be fitted with an orbital solution with 
a period of 711 days and an eccentricity of 0.4. The lower 
panel shows the residual velocity. The rms scatter of the 
residuals is partly due to intrinsic variability on short time 
scales, as commonly found for giant stars. 



function for HD 11977 in a more rigorous way. The 
subsolar metallicity of this star ([Fe/H] = —0.21 ± 0.1) 
lead us to derive a mass of 1.91 ± 0.21 Mq, and to 
exclude masses above 2.5 Mq. A star of such mass and 
metallicity is not expected to have suffered significant 
mass loss during its ascent of the RGB. Therefore, the 
progenitor of HD 11977 was probably an A-type main 
sequence star. 

The projected rotational velocity vsmi was 
computed using the cross-correlation method (e.g., 
Benz & Mayor [TOM Setiawan et al. l2TM}l . From the 
projected rotational velocity and the stellar radius we 
calculated the upper limit of the rotation period as 
P ro t/ sin i = 230-270 days. Such a procedure of P rot 
determination is very sensitive on the v sin i precision. 
Nevertheless, the vsnvi value we have obtained for 
HD 11977 is in excellent agreement with the value 
expected for a solar-type G5III star with B — V ~ 0.90, 
■usini = 2.1 kms" 1 (de Medeiros et al. I1996j> . pointing for 
an excellent precision on our wsini determination. 

3. Observations and analysis 

We have monitored the radial velocity of HD 11977 
from October 1999 to February 2002 with the 
echelle spectrograph FEROS (R = A/ A A = 48 000) at the 
1.52-m ESO telescope. We continued the monitoring at 
the 2.2-m MPG/ESO telescope, to which FEROS had 
been moved, from December 2003 to November 2004. 
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FEROS is equipped with two fibres and can be op- 
erated in "object-sky" and "object-calibration" mode 
(Kaufer & Pasquini I1998[l . We used the simultaneous 
calibration technique (Baranne et al. I199b| and the 
cross-correlation method to determine the radial veloc- 
ity of HD 11977. Details of the data reduction and 
RV computation procedures have been described in 
Setiawan et al. (2003 ). The radial velocity measurements 
of HD 11977 are shown in Fig. 1. 



4. Evidence for a substellar companion 

We observed a long-period variation in the radial veloc- 
ity of HD 11977. A Lomb-Scargle periodogram (Scargle 
1982) of the data yields significant power at a frequency of 
0.0014 cycles day -1 (Fig. 2) and a false alarm probability 
(FAP) of ^10 -9 using expressions given in Scargle H1982J1 . 
This low FAP value was confirmed using the boot- 
strap randomisation technique (Kiirster et al. I1997|) and 
200 000 random shuffles of the data which yielded a 
FAP < 5 x 10~ 6 , i. e. no power in the random data greater 
than the data power. 
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Fig. 2. The Lomb-Scargle periodogram of the radial ve- 
locity measurements. The periodogram shows a significant 
power at a frequency of 0.0014 cycles day -1 . 



It is unlikely that the RV variation of HD 11977 is due 
to rotational modulation. The upper limit of the rotation 
period (P ro t/ sini < 300 days) is obviously much lower 
than the observed period in the RV variability. The anal- 
ysis of the Ca II K emission lines (Fig. 3) did not show 
any significant stellar activity (e.g., due to surface inho- 
mogeneities). We measured the relative flux K2V/K2R. 
This method has been applied on several red giants to 
detect rotational modulation (Setiawan et al. 2005). In 
addition, we measured the bisector velocity span (e.g. 
Queloz et al. I2001|) and found no dependence between the 
asymmetry of the spectral line profile and the measured 
radial velocities (Fig. 4) . 
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Fig. 3. The upper panel shows the region around the 
Ca II K emission lines. The relative flux K2V/K2R (lower 
panel) shows no long-period variation. 



Assuming a semi-amplitude of 105 ms _1 and a rota- 
tional velocity of 2.4 kms -1 , a simple starspot model re- 
sults in a spot-filling factor of 4% of the stellar hemisphere 
(Hatzes 2002). For a starspot close to the equator and 
with AT = 1200 K, this translates to a photometric varia- 
tion of ~0.1 mag (Launhardt & Setiawan, in preparation), 
which is of a factor of ~200 larger than the value given in 
Adelman (2001). However, we cannot exclude the possibil- 
ity that the HiPPARCOS photometric measurements were 
done at a time when HD 11977 was in the "quiet" phase 
of its activity cycle. 

The astrometric variation of HD 11977 is less than 
1 mas (Hipparcos). This gives a lower limit for sini <; 0.1. 
If the sini were 0.1, the true rotation velocity would be 
about 24 kms -1 . A K giant rotating this fast would be 
very active. In this case, strong Ca II K emission lines 
should be expected. We did not observe this feature in 
our spectroscopic measurements. 

We therefore exclude rotational modulation due to sur- 
face inhomogeneities as the source of the RV variation, 
and infer the presence of a substellar companion (planet 
or brown dwarf) orbiting HD 11977. The orbital solu- 
tion to these data yielded a period of P = 711 days, a 
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Fig. 4. Bisector velocity span (BVS) of HD 11977 calcu- 
lated from the spectra taken between October 1999 and 
November 2004. The diagram shows no correlation be- 
tween BVS and the measured radial velocities. 

semi-amplitude of K\ = 105 ms -1 and an eccentricity of 
e = 0.4. The residual of 29.1 ms -1 is comparable with the 
expected long-term accuracy of FEROS. The computed 
orbital parameters yield a minimum mass of the compan- 
ion 77i2 sini = 6.54 Mj up (Table 2). 

Table 2. Orbital parameters for HD 11977 B 



p 


711 ±8 


days 


T (JD - 2450000) 


1420 


days 


e 


0.4 ±0.07 




V 


-18.010 


kms^ 1 


Ull 


351.5 ±9.5 


deg 


Ki 


0.105 ±0.008 


kms -1 


o(0 - C) 


29.1 


DIS" 1 


f(m) 


(6.55 ±1.63) x 10~ 8 


M© 


ai sin i 


9.41 x 10 5 


km 


with mi = =1.91M 






77i2 sin i 


6.54 


M. Iup 


£72 


1.93 


AU 



The HiPPARCOS astrometry data yields an upper limit 
for the companion mass of less than 66 Mj up , i.e. be- 
low or near the hydrogen-burning limit. The expectation 
value of sin i = 7r/4 for randomly oriented orbital planes 
can be used to derive a value of m,2- We derived a value of 
mi p» 8.3 Mj U p. Since this value is below the deuterium- 
burning limit (~13 Mj up ), HD 11977 B is one of very few 
planetary companions discovered around intermediate- 
mass giant stars. 



5. Discussion 

To date there are only two substellar companions de- 
tected around intermediate-mass stars, namely HR 7329 
(Lowrance et al.EESJ and HD 13189 (Hatzes et al. 2005). 



HR 7329 A is a main-sequence star of spectral type 
AO of mass M* = 2.9 M Q . HR 7329 B is a brown- 
dwarf companion with 7712 ~ 30-50 Mj up in an orbit of 
-200 AU (Guenther et al- lSIM) . HD 13189 A is a K gi- 
ant/supergiant of mass = 2-7 Mq. The distance to 
this star is only poorly known, thus its stellar mass cannot 
be determined accurately. This leads to large uncertainty 
for the mass of the substellar companion. The compan- 
ion's minimum mass was estimated between 8-20 Mj up . 
Due to this uncertainty HD 13189 B can be either a planet 
or brown dwarf. 
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Fig. 5. Mass distribution of stars with known planetary 
companions (as of April 2005, including HD 11977). We 
exclude HD 13189 from this histogram due to the large 
uncertainty of the stellar mass. 

HD 11977 is the first intermediate-mass star with ac- 
curate mass around which a planetary companion has 
been discovered (Fig. 5). Assuming negligible mass loss 
during the RGB phase, the main-sequence progenitor of 
HD 11977 is probably an A star with an stellar mass 
of approximately 2 M (Schmidt-Kaler 1982). HD 11977 
has subsolar metallicity ([Fe/H]=— 0.21). Low metal abun- 
dance has been also found in other giants hosting planets 
(Sato et al. 120031 da Silva et al. in preparation). This 
point will be investigated in more detail in our next study 
of stellar abundances. 

Our discovery gives a crucial support to search for 
extrasolar planets around early type main-sequence stars 
and stars in late evolutionary stages. Detections of plan- 
etary companions to stars outside the solar-mass regime 
can provide valuable constraints for theories of planet for- 
mation. 
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